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I.  INTRODUCTION 


The  first  section  of  this  report  deals  with  improvements  to  the  HITRAN  data  base 
that  have  been  implemented  very  recently,  i.e.,  since  the  1986  version  used  for  SHARC 
2.0  was  published.  The  new  O3  data  was  obtained  prior  to  its  publication  in  November 
1992  (Ref.  1).  The  new  1991/92  HITRAN  data  base  includes  many  transitions  that  will 
affect  calculations  for  nonthermal  vibrational  populations.  Several  band  origin  frequencies 
are  also  improved.  In  particular,  the  9.6-10  jim  region  has  been  revised  to  include  many 
new  bands  and  hot  bands.  The  effects  of  the  improvements  to  the  O3  spectroscopic  data 
base  on  SHARC  calculations  of  O3  radiance  have  not  yet  been  quantified;  however,  the 
edges  of  the  fundamental  bands  will  be  affected  the  most,  since  anharmonic  oscillator 
effects  shift  the  band  centers  of  upper  level  transitions  slightly  away  from  those  of  the 
fundamental  bands.  For  V3  >  3,  Av  =  1  emissions,  we  have  tabulated  Einstein  A 
coefficients  derived  by  Rawlins  et  ai.  (Ref.  2)  using  an  anharmonic  oscillator  approach, 
since  there  is  no  information  regarding  those  transitions  in  HITRAN. 

There  is  a  limited  amount  of  experimental  data  pertaining  to  the  mechanisms  and 
rate  constants  for  the  collisional  excitation  and  relaxation  of  ozone  in  the  upper  atmosphere. 
Many  of  the  excitation/relaxation  rate  coefficients  used  in  SHARC  2.0  for  O3  are  based  on 
a  review  article  by  W.  Rawlins  (1985)  (Ref.  3)  and  references  cited  therein.  Rawlins  has 
continued  to  be  active  in  this  subject  area,  with  two  papers  on  the  subject  of  vibrationally 
excited  O3  produced  in  the  three-body  recombination 

O  +  O2  +  M— »C>3  +  M  (1) 

published  in  1987  (Refs.  2,4).  He  has  provided  us  with  helpful  insights  into  the  O3 
scheme  based  on  his  interpretations  of  recent  SPIRIT  I  data  (Ref.  5).  We  have  also 
compared  Rawlins's  proposed  rate  constants  for  reactive  collision  processes  to  those  listed 
in  the  NASA  panel  recommendations  of  1990  (Ref.  6)  and  the  NIST  kinetics  data  base 
(1990)  (Ref.  7). 

Temperature  dependences  of  energy  transfer  processes  involving  O3  are  difficult  to 
study  due  to  the  thermal  decomposition  of  O3.  For  reactive  collision  processes,  the 
temperatures  studied  have  been  less  than  approximately  360  K  (a  temperature  range  more 
relevant  to  the  stratosphere  than  the  upper  atmosphere).  Due  to  the  limited  range  of 
information,  a  T1^2  temperature  dependence  was  used  by  Rawlins  [and  by  Spectral 
Sciences  Incorporated  (SSI),  in  developing  SHARC  2.0]  for  essentially  all  of  the  O3 
energy  transfer  processes  (except  those  with  rate  constants  =  10-11  cm^  molecule*1  s**at 
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room  temperature)  corresponding  to  a  reaction  probability  based  only  on  the  number  of 
collisions. 

After  reviewing  the  rate  expressions  used  in  SHARC  2.0  for  collisional  relaxation 
of  O3,  we  are  still  confused  by  the  apparent  extrapolation  of  a  rate  constant  for  the 
intermode  energy  transfer  from  (00 1 ),( 1 00)- (0 1 0)  to  many  other  vibration-to-translaticn 
(V-T)  processes,  including  several  where  only  a  quantum  of  energy  from  the  V3  mode  was 
relaxed.  The  references  provided  by  SSI  for  this  assumption  refer  to  extrapolations  based 
on  the  harmonic  oscillator  approximation  in  the  SSH  treatment  of  vibrational  relaxation 
rates.  However,  the  extrapolations  made  within  the  SHARC  2.0  O3  scheme  are  not  within 
the  scope  of  the  theory.  In  order  to  remedy  the  current  SHARC  2.0  scheme,  it  appears  that 
first  (100),  (00l)-(000)  V-T  processes  should  be  put  back  into  the  list,  with  estimated  rate 
constants.  The  SHARC  2.0  treatment  of  this  aspect  of  the  O3  scheme  is  very  different  from 
that  of  Rawlins  (Ref.  3),  who  used  the  standard  treatment  of  multiplying  the  (001)-(000) 
rate  constant  by  v  for  upper  level  relaxation  processes  to  v-1. 

One  additional  mechanism,  involving  the  excitation  of  03(102)  and  (101),  seems  to 
be  necessary  in  SHARC  in  order  to  account  for  the  SPIRE  rocket  observations  and 
laboratory  observations  of  Rawlins  et  al.  (Refs.  2,4,5).  A  change  that  would  make  the 
distinction  between  O2  and  N2  as  collision  partners  for  the  O3  (100,001)-(010)  is 
suggested  in  Section  III  of  this  report.  After  reassigning  rate  constants  to  the  V-T  one- 
quantum  energy  transfer  processes  of  the  lowest  vibrational  levels,  the  rest  of  the  O3 
scheme  involving  higher  vibrational  levels  up  to  V3  =  9  can  be  dealt  with. 

A  listing  of  the  current  (SHARC  2.0)  O3  scheme  is  provided  in  Appendix  A.  The 
listing  was  provided  by  SSI  and  differs  from  the  module  actually  in  the  code  only  by  the 
inclusion  of  references.  This  scheme  represents  one  possible  way  of  handling  excitation 
and  relaxation  mechanisms  within  the  complex  array  of  O3  vibrational  energy  levels.  We 
suggest  some  changes  that  may  or  may  not  be  compatible  with  the  overall  assumptions 
made  in  the  SSI  scheme.  Other  updates  that  we  suggest  in  this  report  should  be  relatively 
easy  to  incorporate  into  SHARC  3.0. 
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n.  EINSTEIN  A  COEFFICIENTS 


The  most  fundamental  addition  that  can  be  made  to  the  SHARC  O3  File  is  a  revised 
list  of  vibrational  bands  and  their  spontaneous  emission  coefficients.  There  have  been 
many  improvements  and  additions  to  the  HITRAN  data  base  for  O3  since  the  1986  version 
was  published.  The  improvements  were  achieved  by  a  concerted  effort  of  investigators, 
who  spent  a  considerable  amount  of  time  on  improving  the  assignments  and  band  strengths 
for  the  1991  HITRAN  listing.  A  publication  describing  the  1991/92  HITRAN  data  base 
can  be  found  in  the  Journal  of  Quantitative  Spectroscopy  and  Radiative  Transfer  (Ref.  1). 
The  contributors  to  the  O3  section  are  referenced  therein.  HITRAN  91/92  contains  57 
rovibrational  bands  for  the  main  666  isotope  of  O3  and  nine  sets  of  pure  rotational  lines,  in 
contrast  to  HITRAN  1986,  which  contained  only  16  rovibrational  bands  and  one  set  of 
pure  rotational  transitions  within  the  (000)  manifold.  Many  of  the  improvements  to 
HITRAN  are  relevant  to  SHARC.  For  example,  Vo  for  the  (003)-(000)  band  has  been 
reassigned  as  3046.0082  cm'1  rather  than  the  previous  value  of  3041.200  cm*1.  These 
types  of  improvements  will  necessitate  changes  in  the  03STAT.DAT  file  of  SHARC  as 
well  as  in  the  03KIN.DAT  file.  Some  of  the  band  strengths  previously  estimated  in 
SHARC,  such  as  that  of  the  (200)-(100)  transition,  have  now  been  included  in  HITRAN. 
The  Einstein  A  coefficient  for  that  transition  was  previously  given  as  1.18  s*1  in  SHARC 
2.0,  whereas  using  the  new  HITRAN  band  strength,  A  is  calculated  to  be  0.640  s*1.  This 
is  nearly  a  100%  decrease  in  the  probability  for  spontaneous  emission  from  that  band, 
which  occurs  at  8.62  p.m.  The  (003)-(001)  A  value  has  changed  from  1.03  s'1  to 
0.571  s-1.  The  overall  effect  that  the  updates  will  have  on  radiance  calculations  performed 
with  SHARC  has  not  yet  been  quantified,  and  many  of  the  "new"  transitions  occur  from 
vibrational  levels  that  will  only  be  populated  under  very  non-LTE  conditions.  However, 
since  the  code  requires  that  levels  up  to  as  high  as  V3  =  9  be  populated,  all  of  the  additional 
transitions  should  be  included.  A  listing  of  all  the  O3  bands  listed  in  HITRAN  91/92  (Ref. 
1)  is  given  in  Table  1  with  notations  indicating  changes  from  the  1986  listing  and/or  from 
SHARC  2.0. 

Band  origins  and  Einstein  A  coefficients  for  observed  AV3  =  1  transitions  have  been 
reported  by  Rawlins  et  al.  (Ref.  4)  from  experiments  in  which  vibrationally  excited  O3  was 
generated  by  reaction  (1).  The  chemiluminescence  from  product  (>3(v)  was  observed  at 
relatively  high  resolution  and  at  low  rotational  temperatures.  The  band  centers  retrieved 
from  the  best  fits  to  the  experimental  data  are  given  in  Table  IB. 
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Table  1:  Einstein  A  coefficients  for  O3 


A:  From  the  HTTRAN  database  199 1/92  Edition  (Ref.  1) 

vo(ltxn)  vo  (cm*1)  v’-v”  ES  A  (s'1) 


cm-  V/molcculc-cm'^) 


*164 

61.0533 

100  -  001 

2.039E-22 

0.00035 

*144 

69.7394 

110-011 

1.514E-24 

0.00028 

29.31 

341.1529 

001-010 

1.001E-21 

0.0033 

24.86 

402.2062 

100-010 

2.823E-24 

0.000012 

16.04 

623.3852 

011-100 

2.521E-23 

0.0017 

14.68 

681.4119 

021-011 

1.958E-22 

0.314 

14.61 

684.4385 

011-001 

3.596E-21 

0.209 

14.49 

690.3148 

120-110 

1.181E-22 

0.272 

14.43 

693.1246 

110-  100 

2.624E-21 

0.210 

14.32 

698.3415 

020  -  010 

3.879E-20 

0.445 

14.27 

700.9311 

010-000 

5.921E-19 

0.227 

13.26 

754.1779 

110-001 

2.461E-23 

0.0017 

10.75 

929.8447 

012-110 

3.290E-23 

0.096 

10.47 

954.7535 

002-100 

3.2C5E-21 

0.483 

10.11 

988.9772 

111  - 110 

2.102E-21 

9.68 

10.00 

999.5841 

012-011 

5.996E-21 

20.09 

9.92 

1007.6470 

101  -100 

6.246E-20 

10.27 

9.91 

1008.6618 

021  -020 

1.497E-20 

10.40 

9.87 

1015.8068 

002  -  001 

1.680E-19 

20.83 

9.75 

1025.5914 

011  -000 

4.635E-19 

11.17 

9.60 

1042.0840 

001-000 

1.413E-17 

11.-56 

9.45 

1058.7166 

111-011 

6.379E-23 

0.239 

9.36 

1068.7003 

101-001 

5.238E-21 

0  719 

9.20 

1087.3041 

120  -  020 

2.207E-22 

0.178 

9.18 

1089.9162 

210-110 

3.331E-23 

0.146 

9.13 

1095.3308 

110-010 

1.273E-20 

0.349 

9.11 

1098.0179 

200-  100 

3.299E-21 

0.642 

9.07 

1103.1372 

100  -  000 

5.403E-19 

0.498 

8.62 

1159.0712 

200  -  001 

8.340E-21 

1.34 

8.62 

1159.6556 

210-011 

1.652E-22 

0.742 

7.37 

1356.9597 

002-010 

2.064E-23 

0.00087 

7.15 

1399.2726 

020  -  000 

4.499E-22 

0.0007 

7.09 

1409.8532 

101-010 

9.927E-22 

0.045 

5.86 

1707.0034 

021-010 

3.421E-21 

0.227 

5.79 

1726.5225 

011-000 

5.357E-20 

0.120 

5.60 

1785.6456 

120-010 

1.372E-21 

0.100 

5.57 

1796.2619 

no-ooo 

2.374E-20 

0.058 

5.15 

1942.9509 

003-  100 

1.379E-21 

rs  - r 

U.OJ  / 

5.05 

1980.5651 

102-100 

3.838E-22 

0.242 

4.99 

2004.0042 

003-001 

1.189E-21 

0.571 

4.94 

2025.1755 

012-010 

3.156E-21 

0.295 

4.90 

2041.6184 

102-001 

1.562E-20 

7.78 

4.86 

2057.8908 

002  -  000 

1.120E-19 

0.358 

4.80 

2083.2724 

201  -  100 

1.049E-20 

7.32 

4.80 

2084.3080 

111-010 

4.083E-20 

4.04 

4.74 

2110.7843 

101-000 

1.236E-18 

4.15 

Table  1A  (com.): 

vo()im)  vo(cm_1)  v‘ -v"  X  S  A(s_1) 


4.66 

2144.3257 

201  -001 

4.579E-22 

0.252 

4.58 

2185.2470 

210-010 

7.633E-22 

0.083 

4.54 

2201.1552 

200  -  000 

3.371E-20 

0.123 

4.15 

2407.9345 

021-000 

4.188E-22 

0.0018 

4.02 

2486.5766 

120  -  000 

3.842E-22 

0.0018 

3.67 

2726.1066 

012-000 

3.301E-21 

0.0185 

3.59 

2785.2391 

111-000 

2.425E-20 

0.142 

3.46 

2886.1781 

210-000 

1.101E-21 

0.0069 

3.28 

3046.0882 

003  -  000 

1.312E-19 

0.918 

3.24 

3083.7024 

102  -  000 

1.265E-20 

0.091 

3.14 

3186.4097 

201  -000 

9.300E-21 

0.071 

011-110 

1.447E-24 

001-100 

1.287E-22 

000-000 

4.470E-19 

010-010 

1.449E-20 

001-001 

2.564E-21 

100-100 

1.871E-21 

020-020 

4.582E-22 

011-011 

8.359E-23 

110-110 

5.564E-23 

a)*outside  the  wavelength  region  considered  in  SITARC  2.0 


Table  IB:  From  Rawlins  and  Armstrong  (Ref.  2)  anharmonic  oscillator  treatment. 
vo(jim)  vq  (cm-1)  v’ -v"  A(s'1) 


9.60 

1042 

001-000 

11.(2) 

9.84 

1016 

002-001 

20.7 

9.92 

1008 

101-100 

6.0** 

10.13 

987 

003-002 

28.5 

10.30 

971 

102-101 

17.5 

10.50 

952 

004-003 

34.4 

10.65 

939 

103-102 

24.0 

10.92 

916 

005-004 

38.5 

11.06 

904 

104-103 

28.5 

006-005 

40.0 

105-104 

31.5 

007-006 

40.9 

106-105 

33.0 

b)  italics  indicate  transition  listed  in  HTTRAN 

c)  **  A  significantly  different  from  HTTRAN 


The  five  quantum  level  of  V3  was  the  highest  observed  in  Rawlins's  experiments:  however, 
higher  levels  may  not  have  been  observed  due  to  weaker  intensities  or  resolution 
limitations.  The  authors  discuss  the  possibility  of  detecting  emissions  from  levels  of  V3 
greater  than  5  (which  emit  between  approximately  1 1  and  12  pm)  in  the  upper  atmosphere 
where  number  densities  are  low.  Rawlins  reports  that  recently  analyzed  SPIRIT  I  data 
(Ref.  5),  which  are  of  higher  resolution  than  earlier  field  measurements  (Ref.  8). 
correspond  well  to  the  experimental  data,  with  approximately  the  same  transitions  observed 
(Ref.  2).  For  the  higher  V3  transitions,  the  Einstein  A  coefficients  in  Table  IB  were 
obtained  from  functions  plotted  in  Reference  4,  which  take  into  account  the  anharmonicity 
of  the  oscillator.  For  V3  >  3,  the  A  values  are  reported  to  be  upper  bounds.  Rawlins's 
anharmonic  oscillator  estimation  of  the  Einstein  A  coefficients  is  a  better  one  than  a  simple 
harmonic  oscillator  treatment.  It  is  unclear  where  the  corresponding  coefficients  in 
SHARC  2.0  came  from.  If  they  were  calculated  from  an  anharmonic  oscillator  function 
based  on  lower  transitions  listed  in  the  HITRAN  data  base,  they  are  probably  satisfactory, 
but  they  would  need  to  be  updated  using  HITRAN  91/92  for  the  basis  of  the  calculation.  It 
is  certain  that  they  were  not  taken  directly  from  Rawlins's  original  review  (Refs.  9,10). 

Finally,  two  bands  not  listed  at  all  in  the  HITRAN  data  base  are  included  in 
SHARC  2.0.  In  the  code  a  rate  constant  for  the  hypothetical  transition  (1 1 1)-(001)  is 
equated  to  the  (110)-(000)  value,  and  a  second  hypothetical  transition  (111)-(100)  is 
assigned  a  rate  constant  from  the  (011)-(000)  band.  We  have  not  found  any  published 
observations  or  measurements  for  these  two  transitions  from  the  (111)  level.  There  is  no 
obvious  reason  for  creating  them,  since  anothe-  transition,  (1 1 1)-(010),  has  a  relatively 
strong  Einstein  A  coefficient,  providing  at  least  one  relaxation  channel  from  that  level.  On 
the  other  hand,  several  bands  not  included  in  SHARC  2.0  are  listed  in  Table  1.  For  a  more 
complete  spectral  picture  of  the  9-10  pm  region,  these  additional  bands  in  particular  should 
be  included  in  SHARC. 
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III.  COLLISi'ONAL  EXCITATION  OF  03 


A. 


EXCITATION  VIA  O  +  02+  M  -»  03(v)  +  M 


The  NIST  Chemical  Kinetics  Database  (1990)  (Ref.  7)  has  37  references  to  kinetic 
data  for  this  reaction  over  the  temperature  range  of  approximately  200-360  K  (see  Figure 
1).  The  rate  constant  exhibits  a  steep  negative  temperature  dependence.  The  rate 
coefficients  reported  in  such  data  bases  as  NIST  or  the  NASA  panel  evaluation  are  for  the 
total  thermal  reaction.  Presumably  the  product  distribution  is  over  all  03(v)  states  formed, 
and  the  rate  constant,  therefore,  refers  to  a  summation  over  all  product  states. 

Early  laboratory  experiments  (Ref.  11)  indicated  that  the  excitation  of  03(v)  via 
process  (1)  was  predominantly  in  the  stretching  modes.  In  the  most  recent  experiments 
conducted  by  Rawlins  et  al.  (Ref.  2),  a  further  attempt  was  made  to  sort  out  the  kinetics  as 
a  function  of  the  final  state  distribution.  The  approach  taken  by  Rawlins  (Refs.  2,3)  was  to 
use  a  "prior"  distribution  of  states  according  to  surprisai  theory  to  calculate  the  population 
probability  factors. 


p°(v)  = 


(1-Ev/E)3/2 

X  (1-Ev/EF- 


V 


(2) 


given  the  total  energy  E  available  to  the  system.  The  experimental  work  by  Rawlins 
(1987)  appears  to  substantiate  this  treatment,  used  first  in  his  1985  review  and  later  in 
SHARC  2.0.  The  final  conclusion  (Refs.  2-4)  is  that  the  energy  distribution  of  the  product 
03(v)  formed  in  process  (1)  is  rapidly  randomized  throughout  the  asymmetric  stretching 
mode  03(v3). 

To  obtain  the  appropriate  rate  constants  for  each  of  the  03(v3)  products,  Rawlins 
(Ref.  3)  used  the  rate  expression  for  the  total  reaction  of 


O  +  O2  +  N2  — »  03(Zv3)  +  N2 


(3) 


from  the  work  of  Huie  et  al.  (Ref.  12), 


ktotal  =  1  •  1  x  10'34  e5  i0/T  cm6  molecule'2  s* 1  (4) 
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log<k>  vs  iOOO/T 


821.1  H/LEU  Q  6  71HIP/TR0 
71STU/NIK  A  8  72TCH 
7  2  HU  I /HER  JklO  73BAL/LAR 
73BEU/J0H  Oi2  74R0S/TRA 
74 SHE  AX4  73HIP/SCH 

76H0G/BUR  Vlt  77ARE/SAH 
79ARN/COH  ®18  79ARN/COM 
79ARN/CQM  ■  20  SOKLA/flND 
80KLA/ANO  *22  80KLA/AND 
80KLA/AHD  Q24  80LAL/UER 
80SUG/X  SH  A26  82LIH/LEU 
82LI N/LEU  A28  82LIN/LEU 
84CRO/TRO  <>31  87DEM/GOL 
71STU/NIK  A33  71STU/NIK 
73BAL/LAR  w35  73BAU/LAR 
73BEU/G0H  ®37  73BEU/OOH 
73BEU/JQH  ■  3  9  73BEU/JOH 
89RTK/BAU  *41  09ATK/BAU 
85BOR/COB 


Figure  1:  Rate  constant  data  for  O  +  O2  +  M  ->  O3  +  M  from  the  NIST  Data  Base,  for 
various  collision  partners,  M.  The  data  set  labeled  31  corresponds  to  the  NASA  panel 
recommendation  (unchanged  from  1987  to  1990  versions).  The  data  set  labeled  9  (Huie 
et  al.)  was  depicted  for  M  =  Ar.  Their  results  for  M  =  O2  lie  in  the  middle  of  the  O2,  N2, 
air  data. 
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and  multiplied  it  by  the  relative  population  of  each  of  the  excited  vibrational  levels  obtained 
from  expression  (2).  At  300  K,  the  rate  constant  calculated  from  the  expression  reported 
by  Huie  et  al.  for  M  =  N2  is  essentially  the  same  as  the  value  calculated  from  the  expression 
recommended  by  the  NASA  panel  for  M  =  air, 

ktotal  =  6  x  10'34  (T/300)*2-3  cm6  molecule-2  s'1  (5) 

From  Eqs.  (5)  and  (6),  the  300  K  rate  constant  is  6  x  10'34  cm6  molecule-2  s‘l.  In  the 
Arrhenius  plot  from  the  NIST  data  base,  given  in  Figure  1,  the  expression  of  Huie  et  al.  is 
linear,  while  that  of  the  NASA  panel  is  nonlinear.  Both  sets  of  data  fall  within  the 
"average"  data  for  the  reaction  with  M  =  air,  O2,  or  N2.  In  Figure  2,  rate  constants 
obtained  from  expressions  (4)  and  (5)  in  addition  to  a  linear  fit  through  the  data 
recommended  by  the  NASA  panel. 


ktotal  =  9.27  x  10'35  e562^cm6  molecule-2  s'1. 


are  plotted  as  a  function  of  T.  Unfortunately,  none  of  these  expressions  was  intended  to  be 
extrapolated  to  temperatures  above  approximately  360  K.  From  Figure  2,  it  is  obvious  that 
the  two  different  functional  forms  (linear  and  nonlinear  in  terms  of  log(k)  and  1000/T)  start 
to  deviate  sharply  as  the  temperature  is  increased. 

The  conclusion  from  this  analysis  is  that  for  atmospheric  temperatures  less  than 
approximately  360  K,  either  the  single  data  set  of  Huie  et  al.  or  a  selective  average  of  data 
reported  for  M  =  air  are  nearly  equivalent.  However,  as  temperatures  approach  400  K  and 
above,  a  large  difference  in  rate  constants  occurs  depending  on  the  choice  of  the  function 
for  the  temperature  dependence.  This  should  be  considered  more  carefully  if,  for  example, 
auroral  temperatures  in  the  upper  atmosphere  are  substituted  into  expression  (4)  in 
SHARC. 

B.  ADDITIONAL  EXCITATION  MECHANISMS 

From  data  obtained  in  microwave  discharge  experiments  (Refs.  2,4)  and  from 
SPIRIT  I  (Ref.  5),  Rawlins  found  evidence  for  O3  (101)  and  (102)  transitions  involving 
one  quantum  of  V3.  In  the  laboratory  work,  evidence  for  the  population  of  the  higher 
stretch  combination  levels  (103)  and  (104)  was  also  found.  A  resonant  V-V  energy 
exchange  with  O2  (v  =  2)  has  been  proposed  to  account  for  the  population  of  (102),  which 
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Rate  Constant  (cm0  molecule* 


7.00e-34 

6.00e-34 

S.00e-34 

4.00e-34 

3.00e-34 

2.00e-34 

1 .00e-34 

O.OOe+O 


600  800 
Temperature  (K) 
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Q  Huie  et  al.,  ( M  =  N2),  Linear  Arrhenius  function  (4) 

•  NASA  Panel  Evaluation,  non-linear  function  (5) 

▲  NASA  Panel  Evaluation,  linear  Arrhenius  function  (6) 


Figure  2:  Temperature  dependence  of  rate  constants  for 
0  +  02  +  M-»03  +  M 


can  be  followed  by  relaxation  to  (101).  The  (101)  level  can  then  radiate,  or  undergo  an 
energy  transfer  process.  The  possible  sequence  of  reactions  is 

02(v  =  2)  +  03(000)  ->  02(v  =  0)  +  03(102)  ( 7  ) 

03(102)  +  M  ->  03(101)  +  M  (  8  ) 

Rawlins  et  al.  also  report  that  it  is  likely  that  03(003)  and  subsequently  03(002)  are 
enhanced  to  some  extent  by  the  deactivation  of  O3(102).  With  respect  to  SHARC 
calculations,  the  inclusion  of  steps  (7)  and  (8)  will  not  have  much  effect  on  the  overall 
radiance,  but  the  fine  structure  of  calculated  spectra  will  be  affected. 

Rawlins  suggests  that  population  of  the  (103)  and  (104)  states  in  the  laboratory 
experiments  occurred  by  rapid,  near  resonant  V-T  transfer  from  03(004)  and  (005)  formed 
in  reaction  (1);  however,  this  thermalization  may  have  been  specific  to  the  laboratory 
experiments.  Two  other  mechanisms  described  in  Reference  2  for  population  of 
combination  levels  by  electronically  excited  02,  or  by  energetic  electrons,  are  also  more 
relevant  to  the  discharge  experiments  than  to  atmospheric  conditions. 
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IV.  COLLISIONAL  RELAXATION  PROCESSES  FOR  03(v) 

Before  discussing  some  of  the  individual  rate  constants  for  collisional  relaxation  of 
03(v)  by  O2,  N2,  and  O,  some  assumptions  made  in  the  overall  scheme  used  in  SHARC 
2.0  should  be  examined  The  scheme  does  not  include  the  processes 


and 


03(001)  +  M'  =  03(000)  +  M' 
03(100)  +  M’  =  03(000)  +  M* 


M'  =  02,  N2  or  O 
M'  =  O2,  N2  or  O 


(9) 

(10) 


While  the  Vi  and  v3  fundamental  bands  are  weaker  than  the  V2  band,  this  doesn't  mean 
that  Processes  (9)  and  (10)  may  be  eliminated.  [This  aspect  of  the  SHARC  2.0  scheme 
was  not  based  on  Rawlins's  review,  since  (9)  is  included  in  his  scheme.]  If  one  assumes 
that  the  (100)  and  (001)  levels  rapidly  equilibrate,  then  possibly  the  (001)  V-T  channel, 
given  in  Eq.  (9),  can  be  used  to  represent  both  pathways  (9)  and  (10).  However,  an 
inconsistency  exists  in  the  SHARC  2.0  scheme  because  processes  are  itemized  fcr  both 
(100)  and  (001)  transfer  to  (010).  This  wouldn't  be  necessary  if  all  the  (100)  energy 
relaxes  through  the  (001)  level.  In  the  actual  experiments  done  to  date,  it  has  not  been 
possible  to  distinguish  between  relaxation  from  (001)  and  (100). 

Even  though  very  little  data  exists  regarding  the  direct  V-T  relaxation  from  (001, 
100),  the  rate  constant  has  been  measured  with  03  as  the  collision  partner  (Refs.  13,14). 
Furthermore,  the  upper  limit  for  the  rate  constant  for  processes  (9)  and  (10),  where  M  = 
O2  and  N2,  is  reported  to  be  approximately  1/10  that  of  the  V-T  relaxation  from  O3(010) 
(Ref.  14).  These  processes  are  therefore  not  necessarily  negligible. 

Keeping  the  omission  of  processes  (9)  and  (10)  in  mind,  the  most  confusing  pan  of 
the  SHARC  2.0  scheme  is  that  all  03(00v)-(00v-l)  V-T  rate  constants  are  set  equal  to 
multiples  of  the  03(001)-(010)  rate  constant  For  example,  the  rate  constant  for 


03  (005)  +  M  =  03(004)  +  M 


is  set  equal  to  5  times  the  rate  constant  for 


k  =  1.45  x  10’ 14  TV2 


(11) 


03  (001)  or  (100)  +  M  =  O3(010)  +  M  k  =  2.89  x 


(12) 


Normally,  the  fundamental  process  involving  relaxation  of  one  quantum  of  energy  from  the 
same  vibrational  mode  would  be  used  as  the  basis  for  the  extrapolation.  The  references 
given  by  SSI  in  Appendix  A,  SHARC  2.0,  do  not  explain  the  assumptions  made  here. 
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Since  processes  (11)  and  (12)  involve  transfer  or  exchange  of  different  types  of  vibrational 
quanta,  the  extrapolation  of  the  rate  constant  for  (1 1)  to  that  for  (12)  seems  unusual.  In 
SHARC  2.0  the  assumption  that  the  rate  constant  for  relaxation  of  one  quantum  of  V3 
energy  should  correspond  to  that  for  process  (12)  was  used  for  10  more  processes  in 
which  a  quantum  of  energy  from  one  mode  of  a  combination  level  was  relaxed  down  to  a 
fundamental  level.  Smce  the  processes  listed  did  not  involve  energy  exchange  between  the 
stretch  and  bend  modes,  it  again  seems  unusual  that  the  rate  constant  for  process  (12)  was 
extrapolated  to  so  many  different  pathways  in  the  SHARC  2.0  scheme. 

The  use  of  only  the  (OOv)  part  of  the  ladder  of  vibrational  levels  in  SHARC  2.0 
appears  to  be  an  artificial  way  of  relaxing  the  upper  levels  to  the  fundamentals  without 
including  the  many  intermode  energy  transfer  pathways  available  between  the  stretch  and 
the  bend  modes.  This  is  somewhat  inconsistent  with  the  inclusion  of  energy  transfer 
processes  from  both  (100)  and  (001)  to  (010)  at  the  fundamental  level.  If  the 
approximation  of  only  using  part  of  the  ladder  is  a  good  one,  then  the  approximation  that 
the  conversion  of  (100)  to  (001)  is  rapid  and  essentially  complete  in  one  direction  should 
also  be  maintained. 

The  potential  errors  associated  with  making  inconsistent  assumptions,  such  as  those 
described  above,  are  compounded  when  M  =  O.  Since  very  little  experimental  data  are 
available  for  O  as  the  collision  partner  for  O3  relaxation,  most  of  the  rate  constants 
associated  with  O  in  SHARC  2.0  are  approximations  obtained  by  multiplying  the  rate 
constant  for  a  comparable  process  involving  N2  or  O2  by  100.  Specifics  of  the  processes 
where  M  =  O  are  discussed  in  Section  HIB. 

A .  COLLISIONAL  RELAXATION  OF  03(v)  WITH  M  =  N2  OR  02 

Recent  infrared  double  resonance  studies  by  Menard-Bourcin,  Menard,  and 
Doyennette  (Refs.  13-15)  (1989-1990)  have  provided  measurements  specific  to  N2,  02, 
At,  and  O3  as  collision  partners  in  the  relaxation  of  vibrationally  excited  C>3(100, 001)  and 
03(010).  The  error  bars  associated  with  their  results  are  significantly  improved  over  those 
of  earlier  measurements,  particularly  those  of  West,  Weston,  and  Flynn  (1976)  (Ref.  16). 
One  significant  feature  of  the  experiments  performed  by  Menard-Bourcin  et  al.  (Ref.  14)  is 
that  measurements  were  actually  made  with  N2.  While  O2  and  Ar  apparently  behave 
similarly  in  most  interactions  with  03,  the  intermode  transfer  between  03(100,001)  and 
O3(010)  occurs  -  40%  faster  with  M  =  N2  than  with  M  =  O2.  The  six  relevant  processes 
for  which  Menard-Bourcin  et  al.  (Ref.  14)  have  provided  new  data  are  as  follows: 
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03(100,001)  +  M  =  03(010)  +  M  M  =  O2,  k  =  2.9  x  10-14  cm3  molecule-1  s-1  (13a) 

M  =  N2,  k  =  4.0  x  10-14  cm3  molecule-1  s-1  (13b  ) 


03(010)  +  M  =  03(000)  +  M  M  =  O2,  k  =  3.4  x  10' 14  cm3  molecule'1  s'1  (14a) 

M  =  N2,  k  =  3.1  x  10'14  cm3  molecule-1  s-1  (14b) 

03(100,001)  +  M  =  03(000)  +  M  M  =  O2,  k  <  3.1  x  10-15  cm3  molecule-1  s-5  (15a) 

M  =  N2,  k  <  3.1  x  10-15  cm3  molecule-1  s-1  (15b) 

Several  other  papers  describing  work  in  this  area  were  previously  published.  Joens, 
Burkholder,  and  Bair  (Ref.  17)  used  a  somewhat  indirect  method  for  extracting  rate 
constants  from  vibrational  population  distributions  from  vibrational  spectra  of  O3.  which 
may  have  resulted  in  their  anomalous  (100,001)  to  (010)  rate  constant.  West,  Weston,  and 
Flynn  (Ref.  16)  have  assigned  a  large  uncertainty  to  that  rate  constant.  Beyond  that,  the 
choice  of  which  of  the  rate  constants  with  M  =  O2  are  best  is  somewhat  subjective. 
Averaging  the  values  of  the  three  most  reliable  works-Adler-Golden  and  Steinfeld  (Ref. 
18),  West,  Weston,  and  Flynn  (Ref.  16),  and  Menard-Bourcin  et  al.  (Ref.  14)-for  the  M  = 
O2  processes,  would  make  the  difference  between  M  =  O2  and  M  =  N2  indiscernible. 
Based  on  the  many  detailed  and  reliable  experiments  published  by  Doyennette  et  al.  for 
other  vibrational  relaxation  systems,  we  would  suggest  simply  using  their  rate  constants  for 
the  first  four  processes  listed  above.  An  efficiency  coefficient  to  differentiate  between  M  = 
N2  and  M  =  O2  could  be  used  in  the  SHARC  kinetics  module. 


B .  COLLISIONAL  ENERGY  TRANSFER  FROM  O3  (v)  TO  O  ATOMS 

At  altitudes  up  to  -  120  km,  O2  and  N2  rather  than  O  will  be  the  predominant 
collision  partners  for  C>3(v).  Above  that  altitude  atomic  O  begins  to  predominate,  but  the 
O3  number  density  decreases  very  rapidly.  Therefore,  accurate  assumptions  regarding  rate 
constants  for  relaxation  processes  involving  O  are  somewhat  less  important  than  for  M  = 
O2  or  N2.  There  are,  however,  some  modifications  that  would  make  the  SHARC  scheme 
more  consistent. 

Both  a  reactive  and  a  nonreactive  channel  are  possible  when  C>3(v)  collides  with  O, 
namely. 


03(v)  +  O  — >  02  (16) 

03(v)  +  O ->  03(v-l)+ O  (17) 
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The  kinetic  data  pertaining  to  the  relaxation  of  03(100),  (001),  and  (010)  by  0(3P)  are  less 
conclusive  than  for  V-T  relaxation  by  N2  and  O2.  Two  papers  by  West  et  al.  (Refs. 
16,19),  dating  from  1976  and  1978,  were  referenced  with  respect  to  the  O3  +  0(3P) 
process  in  the  review  article  by  Rawlins  (Ref.  3).  It  appears,  however,  that  the 
conclusions  drawn  by  Rawlins  were  based  solely  on  the  earlier  work.  In  the  latter  paper, 
West  et  al.  showed  that  the  chemical  reaction  channel  (16)  is  a  minor  one  compared  to  the 
nonreactive  relaxation  pathway  (17).  In  SHARC  2.0  and  in  Rawlins’s  review,  it  is 
assumed  that  all  of  the  reactive  channels  originate  from  (OOv)  and  that  all  of  the  V-T 
relaxation  to  the  ground  level  proceeds  from  (010),  i.e., 

03(00v)  +  O  ->  O2  (18) 

03(010)  +  O  O3+  O.  (19) 

In  their  latter  paper.  West  et  al.  (Ref.  19)  reported  that  at  least  70%  of  the  O3  deactivation 
occurred  via  the  nonreactive  channel  (17)  and  that  the  upper  limit  for  the  rate  constant  of 
(16)  is  4.5  x  10*12cm3  molecule*1  s*1,  with  a  total  rate  constant  for  both  (18)  and  (19)  of 
10*11  cm3  molecule*1  s*1.  In  contrast  to  this,  a  rate  constant  for  all  reactive  processes  (18) 
where  v  <  9,  of  1  x  10*11  cm3  molecule*1  s*1,  is  used  in  SHARC  2.0.  This  rate  constant 
appears  to  be  based  on  an  early  approximation  by  West  et  al.  (Ref.  16),  which  was  used  in 
Rawlins's  review  (Ref.  3).  This  interpretation,  used  in  SHARC  2.0,  contradicts  the  later 
work  of  West  et  al.,  since  it  would  mean  that  the  rate  constant  for  (19)  is  zero,  and  they 
found  that  approximately  70%  of  the  reaction  went  by  that  pathway.  Furthermore,  it  is 
likely  that  the  reaction  rate  constant  increases  with  increasing  vibrational  excitation  of  O3. 
and  that  using  a  rate  constant  which  is  fast  or  approaching  gas  kinetic  for  the  lowest 
reactant  level  O3  (v  =  1)  limits  how  much  the  rate  constant  for  higher  energy  levels  can  be 
increased.  In  order  to  define  a  lower  limit  for  the  rate  constant  of  the  reactive  pathway 
(18),  we  considered  the  literature  data  for  the  thermal  O3  +  0(3P)  ->02+02  process. 
The  NASA  panel  review  and  the  NIST  data  base  report  that  the  reaction  proceeds  with  a 
rate  constant  of  approximately  8  x  10*15  cm3  molecule*1  s*1  at  room  temperature  (Refs. 
6,7).  All  of  these  references  indicate  that  the  rate  constant  for  reaction  (18)  with  V3  =  1  is 
actually  smaller  than  the  value  of  1  x  10*11  currently  in  SHARC  2.0  and  favored  by 
Rawlins.  We  recommend  that  a  value  of  approximately  4.5  x  10*12  (i.e.,  the  upper  limit) 
be  used  for  the  rate  constant  of  (18)  when  v  =  1  and  that  the  value  be  scaled  with 
increasing  v.  We  consider  next  the  rate  constant  for  (19).  In  SHARC  2.0,  the  rate 
constant  for  relaxation  of  the  (010)  level  when  M  =  O2  or  N2  was  simply  multiplied  by  100 
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to  obtain  a  value  of  3  x  10'12cnrt^  molecule'1  s'1  at  300  K  for  (19).  Coincidentally,  this 
would  be  consistent  with  our  suggestion  of  reducing  the  rate  constant  for  (18)  so  that  the 
rate  constants  for  (18)  and  (19)  totaled  less  than  1  x  10* 1 1  cm^  molecule'1  s'1.  Therefore, 
the  rate  expression  for  (19)  could  be  left  as  it  is. 

In  the  list  of  references  furnished  to  us  by  SSL  the  rate  constant  for  the  energy 
transfer  from  either  of  the  stretching  modes  to  the  bend, 

03(100)  or  (001)  +  O  =  03(010)  +  O  k  =  9  x  10'12  cm3  molecule*1  s'1  at  300K  (20) 

or  k  =  5.2  x  10* 12  T1/2, 

is  referenced  to  the  1976  paper  by  West  et  al.  The  value  of  approximately  10'11  cm^ 
molecule1  s*1  is  basically  an  assumption  made  in  order  to  substantiate  the  two-step 
relaxation  pathway,  (001,  100)-(010)-(000),  for  which  the  latter  step  is  presumably  rate 
limiting.  A  value  of  3  x  10"12  cm^  molecule*1  s*1  for  the  rate  constant  of  the  second  step, 
process  (19),  is  used  in  SHARC  2.0.  After  examination  of  the  SHARC  2.0  scheme,  it  is 
not  clear  why  the  assumption  about  the  rate  constant  for  (20)  is  multiplied  and  extrapolated 
many  more  times  in  the  scheme  for  processes  that  involve  V-T  relaxation  of  V3  rather  than 
an  intermode  energy  transfer.  Again,  there  is  no  (100)  or  (001)  V-T  transfer  to  (000)  listed 
in  the  SHARC  2.0  scheme  for  M  =  O,  which  would  seem  a  more  logical  basis  for 
extrapolations  to  processes  involving  higher  vibrational  levels.  This  is  analogous  to  the 
extrapolations  for  M  =  O2  and  N2  discussed  at  the  beginning  of  this  section. 

C.  TEMPERATURE  DEPENDENCES  OF  THE  RATE 
CONSTANTS  FOR  NONREACTIVE  COLLISIONAL 
RELAXATION  PROCESSES 

There  are  no  experimental  data  available  for  the  temperature  dependences  of  non¬ 
reactive  collisional  relaxation  processes  of  vibrationally  excited  O3.  For  the  most  pan, 
Rawlins  (Ref.  3)  assumed  a  T ^2  function,  based  on  a  simple  increase  in  the  number  of 
collisions  as  the  temperature  is  increased.  Predicting  different  slopes  for  the  temperature 
dependences  is  difficult  without  further  experimental  data  at  temperatures  higher  than  room 
temperature. 
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APPENDIX  A 


This  appendix  contains  a  listing  of  the  current  O3  scheme  used  in  SHARC  2.0.  It 
was  prepared  by  Spectral  Sciences  Incorporated,,  and  they  have  provided  the  references  at 
the  end- of  the  table. 
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TABLE  A1 .  REACTIONS  AND  RATE  CONSTANTS  FOR  Q^  MODEL 


Reaction 

Rate  Constant* 

Reference 

O  +  02  +  H  -»  O3  (OOO)  +  M 

2 . 79xlO”3SeSi0  *  °/T 

a,b 

O  +  O2  +  H  -*  03(001)  +  H. 

2.35xl0“3SeS10'°/T 

a,b 

O  +  02  +  K  -»  O3  (002 )  +  M 

1.85xlO-3Se510-°/T 

a,b 

O  +  02  +  M  "♦  03  (003)  +  V. 

1.42xlO“3se510‘0/T 

a,b 

O  +.  02!  +  M  ->  O3  (p04 )  +  K 

1.06xlO~3Se510'0/T 

a,b 

O  +  02  +  M  -*  O3  (005 )  +  M 

7 .59xl0“36e510*°/T 

a,b 

~  +  02  +  M  -»  O3  (006)  +  M 

4.73xlO“36e510-°/T 

a,b 

0  +  02  +  M O3  (007 )  +  M 

2.53xlO“36e510-°/T 

a,b 

0  +  02  +  M  -»  O3  (008 )  +  H 

8.80xl0“37eS10*°/T 

a,b 

O  +  02  +  H  -»  O3  (009 )  +  H 

2.2OxlO“37eS10*°/T 

£ 

O3(010)  +  K  03(000)  rf  M 

1 . 73xlO— — “T1/7 

c 

•o3(oi0)  .+  0  «  o3(ooo)  +  0 

1.73x20_13T1/2 

d 

O3  (020)  +  M  -  O3  (010)  +  M 

3.46x10”1£T1/'2 

e 

■  O3  (020)  +  O  =  O3(010)  +  0 

3 . 46x10“ 13T-/7 

e 

O3  (100)  +  H  .=  03(  OlO )  +  M 

2  •  8SxlO“-£T~/2 

c 

03(100)  +  0  -  03(010)  +  0 

5.20x10~13T1/2 

d 

O3(100)  +  M*  =  03.{001)  +  H* 

1.00x10“ 11 

b,c 

03  (200)  +  K  -  P3(100;  +  H 

S.80xl0"1ST1/2 

e 

O3  (200)  +  O  =  O3(100j'+  O 

1.04x10“12T1/2 

e 

O3  (200)  +  K'  »  O3(101)  +  M 

1 . 00x10“ 2-1 

f 

03(001)  +  H  =  O3(010)  +  K 

l.SSxlO"15!1/2 

c 

03(001)  +  0  .=  O3(01Q)  +  0 

.£ .  20x!0“-3T3/2 

d 

03  (002 )  +  H  =  03(001)  •+■  M 

S.80x10“1£T1/2 

e 

03  (002 )  +  0  =  03(001)  +  0 

1.04x10“2-2T1/2 

e 

O3  (003 )  +  M  =  O3(002)  +  H 

S.70x10“15T1^2 

e 

o3  (003 )  +  0  =  03(002)  +  a 

1.S6x10"12T1/2 

e 

03  (004  )  +  K  =  03 (003 )  +  K 

1.16x10“74T2'/2 

e 

03(004)  +  0  =  O3  (003 )  +  0 

2 . 08x10“ 32 T7/2 

e 

O3(00S).+  H  =  O3  (004 )  +  H 

l^SxlO*-2-4?1/2 

e 
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TABLE  A1  .  (CONTINUED) 


Reaction 

O3  (005)  +  O  =  03(004 )  +  O 
O3  (006)  +  M  a  O3  { 005  )  +  H 
03  (006)  +  O  =  03  (005 )  +  O 
O3  (007)  +  M  »  03  (006)  +  M 
03  (007 ):  +  O  =  03  (006 )  +  O 
03  ( 008}  +  M  =  O3  (007 )  +  M 
03  (008)  +  O  «  03  (007 )  +  O 
03(009)  +  H  ~  O3  (008)  +  M 
O3  (009 )  +  O  =  O3 (008)  +  O 
03  (.001)  +  0  ,02  +  O2 

03  (002)  +  O  -»  02  +  O2 
03  ( 003 )  +  O  -*  02  +  02 
O3(004)  +  O  -»  02  +  02 
03  (005)  +  O  -»  o2  +  o2 
03  (006 )  +  O  -»  02  +  02 
03  ( 007 )  +  O  02  +  02 
O3  (008 )  +  O  -*  02  +  02 
03(009)  +  O  -»  02  +  o2 
O3(011)  +  M  »  03(010>  +  K 
O3(011)  +  O  =  03(010)  +  O 
O3(011)  +  H  =  03(001)  +  M 
O3(011)  +  O  =  O3(001)  +  O 
O3(110)  +  M  =  03(010)  +  H 
O3(110)  +  O  =  03(010)  +  O 
O3(110)  +  H  =  O3(100)  +  H 
O3(110)  +  O  =  O3(100)  +  O 
O3(101)  +  H  =  O3(001)  +  H 

o3(ioi)  +  o  =  o3(ooi)  +  o 

O3(101)  +  K  =  O3(100)  +  H 


Rate  Constant* 

Reference 

2.60xl0~12T1/2 

e 

1.74X10-14!1/2 

e 

S.^xlO-12!1/2 

e 

2.03x10-14T1/2  ' 

e 

Z.ttxlQ-12?1/2 

e 

2.32xl0-14T1/2 

e 

4.16x10-12T1/2 

e 

2.61x10-14T1/2 

e 

4.68x10“12T1/2 

e 

l.OOxlO-11 

b,d 

l.OOxlO-11 

b,d 

l.OOxlO-11 

b ,  d 

l.OOxlO-11 

b,d 

l.OOxlO-11 

b,  d 

l.OOxlO-11 

b ,  d 

l.OOxlO-11 

b ,  d 

l.OOxlO-11 

b ,  d 

l.OOxlO-11 

b,d 

2.89xlO-lsT1/2 

5 

5 . 20xl0-11T1/2 

5 

l^xlO-15!'1/2 

5 

1.73x10-12T1/2 

9 

l.SSxlO-15!1/2 

g 

5 •20xl0-12T1/2 

g 

1.73xl0-lsT1/2 

g 

l^xlO-13?1/2 

g 

2.89xlO-lsT1/2 

g 

S.20xl0-13T1/2 

g 

2.89xlO-lsT1/2 

g 
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TABLE  Al.  (CONTINUED) 


Reaction  Rate  Constant*  Reference 


O3{101)  +  0  -  03(100)  + 

0 

S.ZOxlO-^T1/2 

g 

03(101)  +  H*  =  03 (002)  +  M' 

l.OOxlO-11 

f 

o3(iii)  +  a  —  o3(ioi)  + 

H 

1 . 73xlO-33T3/2 

g 

o3(iii)  + -o  —  o3(ioi)  + 

0 

l.TlxlO-13!1/2 

g 

03(111)  +  M  -  03(011)  + 

M 

2.89xio~15t1/2 

g 

03{111)  +  0  =  03 (Oil)  + 

o 

S^OxlO-13!1/2 

5 

03(111)  +  M  •=  03(110)  + 

M 

2.89xl0-15T1/2 

g 

,03(111)  +  O  =  O3(110)  + 

0 

s^oxio-13!-1/2 

g 

03(010)  -♦  O3(000)  +  hi/ 

0.224 

h 

03(001)  -»  03(000)  +  hi/ 

11.00 

h 

03(100)  -*  O3(.000)  +  hi/. 

0.594 

h 

O3(011)  -*  O3(000)  +  hi/ 

0.116 

h 

O3(110)  -*  03(000)  +  hi/ 

0.053 

h 

03  (002 )  -*  O3(000)  +  hi/ 

0.342 

h 

o3(ioij  ->  o3(6oo)  +  hi/ 

3.67 

h 

03 (200)  -»  03(000)  +  hi/ 

0.106 

h 

03(111)  -*  O3{000)  +  hi/ 

0.131 

h 

03  (003)  -»  03(000)  +  hi/ 

0.746 

h 

03  (020)  ->  O3(0.10)  +:  hi/ 

0.446 

h 

03(011)  -*  O3(010)  +  hi/ 

10.39 

h 

O3(110)  -*  O3(010)  +  hi/ 

0.291 

h 

03(111)  -*  O3(010)  +  hi/ 

3.67 

h 

03  ( 002 )  -*  O3(001)  +  hi/ 

20.72 

h 

O3(101)  -»  O3(001)  +  hi/ 

0.437 

h 

03(111)  -*  O3{001)  +  hi/ 

0.053 

h 

03  (003)  -♦  O3(001)  +  hi/ 

1.03 

h 

O3{101)  O3(100)  +  hi/ 

9.68 

h 

03  (200)  -♦  O3{100)  +  hi/ 

1.18 

h 

03(111)  03(100)  +  hi/ 

0.116 

h 
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TABLE  Al.  (CONTINUED) 


Reaction 


Rate  Constant*  Reference 


o3{Hi) 

-»  03  (011) 

+ 

hi/ 

0.593 

h 

03(111) 

-4  03(110) 

+ 

hv 

'll.Ol 

h 

O3  (003 ) 

-+  O3  (002  ) 

+ 

h  v 

30.44 

h 

03(004) 

-4  O3  (003  ) 

+ 

hu 

36.34 

h 

O3(00S) 

-♦  03  (004 ) 

+ 

hu 

43.08 

h 

C3 (006) 

-4  O3  (005  ) 

+ 

h:/ 

48.49 

h 

03 (007 ) 

-4  03  (006) 

+ 

hi/ 

52..  71 

h 

03 (008) 

-4  03  (007  ) 

+ 

hi/ 

'56.12 

h 

03 (009 ) 

-4  03  (008) 

+ 

hv 

58.67 

h 

Rate  constants  are  in  units  of  1/s  for  unimolecular  reactions, 
cm^/molecule-s  for  bimalecular  reactions  and  cm6/molec,ule-s  for 
termolecular  reactions. 

a  R.E.  Huie,  J.T.  Herron,  and  D.D.  Davis,  "Absolute  Rate  Constants 
for  the  Reaction  of  O+O2+M-4O3+M  over  the  Temperature  Range  200-346 
X,"  J.  Phvs .  Chem. ,  76,  2653  (1972). 

k  W.T.  Rawlins,  "Chemistry  of  Vibrationally  Excited  Ozone  in  the 

Upper  Atmosphere,"  J.  Geoohvs.  Res..  SO,  12283  (1985). 

c  S.M.  Adler-Golden  and  J.I.  Steinfeld,  "Vibrational  Energy  Transfer 
in  Ozone  by  Infrared-Ultraviolet  Double  Resonance,"  chem.  Phvs. 
Letters,  76,  479  (1980). 

d  G.A.  West,  R.E.  Weston,  Jr.,  and  G.W.  Flynn,  "Deactivation  of 

Vibrationally  Excited  ozone  by  0{3p)  Atoms,"  Chem.  Phvs.  Lett.,  42 , 

•  488  (1976). 

e  Estimate  based  on  SSH  theory,  kn_,n_3  =  nk^Q. 

f  Estimate  based  on  Mmn-l-lron+l  =  k100-+00l* 

9  Estimate  based  on  klmn-»lm-ln  =  *010-000'  klmn-*l-lmn  =  kl00-000'  or 
klmn-*lmn-l  =  k001-*000  • 

h  L.  S.  Rothman,  et  al. ,  "The  HITRAN  Database:  1986  Edition,"  Aool . 
Oot. .  16,  4058  (1987). 
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